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Abstract: I review the NLO QCD calculations of the branching ratio for B Xs^ in the SM. 
Including the leading electromagnetic corrections, one obtaines BR{B Xs^) — (3.32 ±0.30) x 10""*. 
Confronting theory with the newest data, an updated value for \Vts\ is obtained: \Vts\ = 0.037± 0.007. 
Theoretical progress on the photon energy spectrum is also discussed. The inclusive FCNC semilep- 
tonic decays in the SM are briefly summarized. Furthermore, B ~+ Xs^ is considered in 2HDMs and 
in different SUSY scenarios. QCD corrections are shown to be crucial. 



> 

oo 



ix. 



X 



1. Introduction 

In the Standard model (SM), rare B meson de- 
cays are induced by one-loop diagrams, where 
W bosons and up-type quarks are exchanged. In 
many extensions of the SM, there are additional 
contributions, where the SM particles in the loop 
are replaced by nonstandard ones, like charged 
Higgs bosons, gluinos, charginos etc. Being in- 
duced also at the one loop-level, the new physics 
contributions are not necessarily suppressed rel- 
ative to the SM one. The resulting sensitivity 
for nonstandard effects implies the possibility for 
an indirect observation of new physics, or allows 
to put limits on the masses and coupling param- 
eters of the new particles. A general overview, 
where many B physics observables are investi- 
gated w.r.t. new physics, was given at this con- 
ference by A. Masiero Concerning the new 
physics aspects of rare B decays, I concentrate in 
this article (see sections 3 and 4) on recent cal- 
culations of the branching ratio for the process 
B — > in a general class of two-Higgs-doublet 
models (2HMDs) § | and in supersymmetric 
scenarios [§, |^, ^ , stessing in particular the im- 
portance of leading- (LO) and next-to-leading 
(NLO) QCD corrections. 

However, also in the absence of new physics, 
rare B decays are very important; they can be 
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used for the determination of various CKM ma- 
trix elements, occurring in the SM Lagrangian. 
To extract these parameters from data, it is cru- 
cial that the corresponding decay rates are reli- 
ably calculable. An important class of such de- 
cays are the inclusive rare B decays, like B 
X,7, B ^ XJ+£-, B ^ Xsv9, B ^ X^-f, 
B Xdi'^i^, B Xdvv, which are sensitive 
to the CKM matrix elements |Vts| and \Vtd\, re- 
spectively. In contrast to the corresponding ex- 
clusive channels, these inclusive decay modes are 
theoretically cleaner, in the sense that no specific 
model is needed to describe the final hadronic 
state. Nonperturbative effects in the inclusive 
modes are well under control due to heavy quark 
effective theory. For example, the decay width 
T{B ^s7) is well approximated by the par- 
tonic decay rate T{b Xgj) which can be an- 
alyzed in renormalization group improved per- 
turbation theory. The class of non-perturbative 
effects which scale like I /ml is expected to be 
wen below 10% [0. This numerical statement 
also holds for the non-perturbative contributions 
which scale like 1 /m^ [|[ |j . 

The framework and the NLO theoretical re- 
sults for the branching ratio of the decay B 
Xs7 are discussed in section 2.1; the photon en- 
ergy spectrum and the partially integrated branch- 
ing ratio for this process are reviewed in section 
2.2; an updated value for the CKM matrix ele- 
ment iVtsl, extracted from the most recent mea- 
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surements of BR{B — > Xg^) and the correspond- 
ing calculations, where also the leading electro- 
magnetic corrections are included, is given in sec- 
tion 2.3; the other inclusive rare decays men- 
tioned above, are briefly discussed in section 2.4. 

The exclusive analogues, B K*j, B 
K^*H+l-, B ^ pj, B ^ pi+£- etc., require the 
calculation of form factors. As the QCD sum rule 
calculations and the lattice results for these form 
factors were summarized by V. Braun [0, I do 
not discuss these decays in the following. 

Finally, there is the class of non-leptonic two- 
body decays, like B tttt, B Ktt; the theo- 
retical status of these processes was discussed by 
L. Silvestrini and new CLEO results were 
presented by D. Jaffe 

2. Inclusive rare B meson decays in 
the SM 

2.1 BR{B Xsj) at NLO precision 

Short distance QCD effects enhance the partonic 
decay rate r(6 —f s"f) by more than a factor of 
two. Analytically, these QCD corrections contain 
large logarithms of the form a" {mi,) log™ {nib /M) , 
where M = mt or M — mw and m < n (with 
n = 0,1,2,...). In order to get a reasonable 
prediction for the decay rate, it is evident that 
one has to resum at least the leading-log (LO) 
series (m = n). As the error of the LO re- 
sult ||l^ was dominated by a large renormaliza- 
tion scale dependence at the ±25% level, it be- 
came clear that even the NLO terms of the form 
as{mb) (a"(TOf,) ln"(mb/M)) have to be taken 
into account systematically. 

To achieve the necessary resummations, one 
usually contracts in a first step an effective low- 
energy theory and then resums the large loga- 
rithms by renomalization group techniques. The 
low energy theory is obtained by integrating out 
the heavy particles which in the SM are the top 
quark and the VF-boson. The resulting effective 
Hamiltonian relevant for b sj in the SM and 
many of its extensions reads 



4Gf 
" V2 



\tY.C,{p)0,{p) 



where Oi (/i) are local operators consisting of light 
fields, Ci{p) are the corresponding Wilson coeffi- 
cients, which contain the complete top- and W- 
mass dependence, and At = VthV^^ with Vij be- 
ing the CKM matrix elements. The CKM depen- 
dence globally factorizes, because we work in the 
approximation A„ = 0. 

Retaining only operators up to dimension 6 
and using the equations of motion, one arrives at 
the following basis 

01 = {cLpl'^bLa) {sLalt^CLp) , 

02 = {cLal^bLa) {sLplpiCLp) , 



O7 



Or = 



167r- 

9s 



A 



-s^ a^^ [mb{p)R) '-f bp G% (2.2) 

As the Wilson coefficients of the QCD penguin 
operators O3, Og are small, we do not list them 
here. 

It is by now well known that a consistent 
calculation for 6 — + S7 at LO (or NLO) precision 
requires three steps: 

1) a. matching calculation of the full standard 
model theory with the effective theory at 
the scale p, = pw to order a" (or a^) for 
the Wilson coefficients, where pw denotes 
a scale of order M\y or mt ; 

2) a renormalization group evolution of the 
Wilson coefficients from the matching scale 
Pw down to the low scale pb = 0{mb), us- 
ing the anomalous-dimension matrix to or- 
der al (or al); 

3) a calculation of the matrix elements of the 
operators at the scale p = pb to order a° 
(or a]). 

At NLO precision, all three steps are rather 
involved: The most difficult part in Step 1 is the 
order as matching of the dipole operators O7 and 
Os- Two-loop diagrams, both in the full- and 
in the effective theory have to be worked out. 
This matching calculation was first performed 
by Adel and Yao flj] and then confirmed by 
Greub and Hurth [ |l5[ , using a different method. 
Later, two further recalculations of this result 

The order ai anoma- 



were presented 



16 



(2.1) lous matrix (Step 2) has been worked out by 
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Chetyrkin, Misiak and Miinz |T^. The extrac- 
tion of certain elements in this matrix involves 
the calculation of pole parts of three loop dia- 
grams. Step 3 consists of Bremsstrahlung con- 
tributions and virtual corrections. The Brems- 
strahlung corrections were worked out some time 
ago by Ali and Greub jisj and have been con- 
firmed and extended by Pott [^. An analysis 
of the virtual two loop corrections was presented 
by Greub, Hurth and Wyler 

Combining the NLO calculations of these 3 
steps, leads to the following NLO QCD predic- 
tion for the branching ratio 0: 



sented in 



we then obtained 



.0.01 _|_ O.OOi 0.29\ rv- 
^ 0.08^ 0.27) ^ 



BR{B^Xsi) = (3.57 ±[,12 

(2.3) 

The central value is obtained for = 4.8 GeV, 
jjLw = TUw and the central values of the input pa- 
rameters listed in |2| . The first error is obtained 
by varying /if, in the interval [2.4, 9.6] GeV, the 
second one by varying the matching scale ^iw be- 
tween mw and m* ; the third error reflects the un- 
certainties in the various input parameters. Sim- 
ilar results were also obtained in refs. |l6|. 



We should mention that in the result (2.3) 
also power corrections are included: there are 
1/m^ corrections 0, whose impact on BR{B 
Xs7) is at the 1% level, as well as nonpertur- 
bative contributions from cc intermediate states 
which scale with 1/™^- Detailed investigations 
H, ^ show that these l/m^ corrections enlarge 
the branching ratio by ~ 3%. 

After the QCD analyses, several papers ap- 
peared where different classes of electroweak cor- 
rections lip, 1^ to BR{B Xs7) were 
considered. In [23|, corrections to the Wilson co- 
efficients at the matching scale due to the top 
quark and the neutral Higgs boson were calcu- 
lated and found to be negligible. The analysis 
PH concluded that the most appropriate value of 

to be used for this problem is the fine struc- 
ture constant a"^ = 137.036 instead of the value 
Q^em = 1 30. 3± 2. 3 previously used. In U|, the 
leading logarithmic QED corrections of the form 
a\og{^w I l^'b) (os log(/iiv//^b))" (with resumma- 
tion in n) were given. 

In ref. |^ we updated the result in eq. ( |2.3| ), 
by including the class of QED corrections pre- 



BR{B ^ X,7) 



/o oo I 0.00 I 0.00 I 0.2&\ -, r)-4 

;^o.o^ in in Q Qgzn Q 25 j X iu 

(2.4) 

The bulk of the change with respect to the value 
in eq. (2^) is due to the different value of a^r^ 
used. 

A remark concerning the error due to the 



variation of the low scale Hb h^ the results (2.3) 
and ( ^.4[ ) is in order here: As it will be dis- 
cussed in more detail in section 3, it was real- 
ized in ||] that in multi Higgs doublet models 
the QCD corrections in certain regions of the pa- 
rameter space are much larger than in the SM. 
As a consequence, the dependence on the scale 
Hb of BR{B — > Xs'y) is also larger. Later, Ka- 
gan and Neubert pointed out very explicitly in 
their analysis |2^ that the scale dependences in 
individual contributions to the branching ratio 
in the SM are larger than their combined effect, 
due to accidental cancellations. They suggest 
that one should add the scale uncertainties from 
the individual contributions in quadrature, in or- 
der to get a more reliable estimate of the trun- 
cation error. Their estimate for the /if, depen- 
dence oiBR{B —> Xs7) is ±6.3%, i.e., more than 
twice the naive estimate. The total error in (2.4), 
however, is dominated by parametric uncertain- 
ties and therefore gets increased only marginally 
when using this more conservative estimate of the 
scale uncertainties. 

The measurement of BR{B Xgj) by the 
ALEPH collaboration at LEP |||| 

BR{B X,-f) = (3.11 ± 0.80 ± 0.72) x 10"'' 

(2.5) 

and by the CLEO collaboration at CESR (H 

BR{B Xs-f) = (3.15±0.35±0.32±0.26)xl0-^ 

(2.6) 

are in good agreement with the NLO calculation 



(2.4), where the most important electromagnetic 
corrections are included. 

2.2 Partially integrated branching ratio in 

B^Xs-f 

The photon energy spectrum of the partonic de- 
cay b S7 is a delta function, concentrated at 
~ (mf,/2), when the b-quark decays at rest. This 
delta function gets smeared when considering the 
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inclusive photon energy spectrum from a B me- 
son decay. There is a perturbative contribution 
to this smearing, induced by the Bremsstrahlung 
process b sjg , as well as a non pertur- 

bative one, which is due to the Fermi motion of 
the decaying b quark in the B meson. 

For small photon energies, the 7-spectrum 
from B — > Xgj is completely overshadowed by 
background processes, like b cud^y and b 
uudj. This background falls off very rapidly with 
increasing photon energy, and becomes small for 
E't- >~ 2 GeV [|||. This implies that only the 
partial branching ratio 

BR{B ^ x,7)(£;r") = ^'^^^ 

(2.7) 

can be directly measured, with i?™™ = 0(2) 
GeV. Recently, CLEO was able to reduce -B™™ 
from 2.2 GeV to 2.1 GeV [|7|. To determine 
from such a measurement the full branching ra- 
tio for B — > Xs"f, one has to know from theory 
the fraction R of the B Xgj events with pho- 
ton energies above Based on calculations 
by Ali and Greub of the photon energy spec- 
trum within the Fermi motion model by Altarelli 
et al. H, CLEO used the value R = 85% - 94% 
Iprt in order to determine BR{B Xg^) from 
the measured partial branching ratio. 

A modern way - based on first principles - 
implements the Fermi motion in the framework 
of the heavy-quark expansion. When probing the 
spectrum closer to the endpoint, the OPE breaks 
down, and the leading twist non-perturbative cor- 
rections must be resummed into the B meson 
structure function f{k^) where fc+ is the 

light-cone momentum of the b quark in the B 
meson. The physical spectrum is then obtained 
by the convolution 

^= r rffc+/(fc+)%^K) > (2.8) 

where {dV part / dE^){m'l) is the partonic differen- 
tial rate, written as a function of the "effective 
mass" ml = rrib + k+. The function f{k+) has 
support in the range — 00 < fc+ < A, where A = 
ruB — TTLi, in the infinite mass limit. This implies 
that the addition of the structure function moves 
the partonic endpoint of the spectrum from rnh/2 



to the physical endpoint itib/^. While the shape 
of the function f{k+) is unknown, the first few 
moments An — J dk+ k"f{k^) are known: Aq = 
1, Al =0 and A2 = -Ai/3. The values of A 
and Al are not calculable analytically and have 
to be extraced from experiments or calculated 
on the lattice. A recent analysis gives (A, Ai) = 
(0.39 iO.llGeV, -0.19 TO.lOGeV^) As A^ 
(n > 2) are poorly known, several Ansatze were 
used for f{k+); e.g. Neubert and Kagan Q used 
f{k+) = N{1 - x)''e(i+'')^, with X ^ k+/A. Tak- 
ing into account the constraints from Aq, Ai and 
A2, the independent parameters in this Ansatz 
can be chosen to be rrib and Ai. As shown in 
||2^ , the uncertainty of mt, dominates the error 
of the partial branching ratio. In Fig. |l| 



4. 




Figure 1: Partially integrated branching ratio as 
a function of the energy cutoff E^^^"; Figure taken 
from Neubert and Kagan |p2| . 

the partial branching ratio is shown for the rel- 
evant range of nib as a function of E'™*", keep- 
ing Ai/A^ fixed. For comparison, the data point 
BR{B Xsj) = (2.04 ± 0.47) x 10-"*, obtained 
in the original CLEO analysis with a cutoff at 2.2 
GeV [Q, as well as the new data point BR{B —^ 
Xs7) = (2.97 ± 0.33 ± 0.30 ± 0.21) x 10-^ cor- 
responing to a cutoff at 2.1 GeV are also 
shown in Fig. ^ We would like to stress that 
it would be very welcome, if the cutoff could be 
pushed down to 2.0 GeV, because the theoretical 
uncertainity drastically becomes smaller, as seen 
in Fig. 

We should add that in the endpoint region 
also the perturbative contributions, mentioned 
at the beginning of this section, could become 
problematic in principle, due to the presence of 
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threshold logs, log(l — 2E^/mh). These loga- 
rithms have recently been resummed up to next- 
to-leading logarithmic precision js^ . The au- 
thors of this paper conclude that for the present 
energy cut at 2.1 GeV, these threshold logs do 
not form a dominant sub-series and therefore their 
resummation is not necessary for predicting the 
decay rate. Recently, also the BLM type cor- 
rections of the order alfi^j to the photon energy 
spectrum were calculated Q. The result was 
used to extract a value for A from the average 
(1 — 2Ej/mB), and a value for Ai from ((1 — 
2Ej/mB)'^)- According to this analysis, the CLEO 
data in the region E-^ > 2.1 GeV implies the cen- 
tral value Aa, ~ 390 MeV and A„2^j, ~ 270 MeV 
at order as and a^Po^ respectively. This anlysis 
was somewhat critisized in ref. js^, by pointing 
out that other terms could be larger than the 
BLM terms. 

2.3 \Vts\ form B^X,7 

Instead of making a prediction for BR{B Xg-f), 
one can use the NLO calculation to extract the 
CKM combination l^tbV't* l/lKbl from the mea- 
surements; in turn, one can determine |Vts|, by 
making use of the relatively well known CKM 
matrix elements Vet and Vtb- Using the CLEO 
(U) and ALEPH data O), one obtains |1 



\Vcb\ 



0.93 ± 0.09e.p. ± O.OSt,, . 



(2.9) 



Using \Vtb\ = 0.99±0.15 from the CDF measure- 
ment and \ Vcb\ = 0.0393 ± 0.0028 extracted from 
semileptonic B decays, one obtains 



WtJ = 0.037 ±0.007, 



(2.10) 



where all the errors were added in quadrature. 
This is probably the most direct determination 
of this CKM matrix element, as the measurement 
oi t ^ sW~^ seems to be difficult. With an im- 
proved measurement of BR{B — > Xg^) and Vtb, 
one expects to reduce the present error on \Vts\ 
by a factor of 2 or even more. 

2.4 B ^ X,a, B ^ {Xs,Xd)l+l- and B ^ 
{Xg,Xd)vv in the SM 

The decay B X^j can be treated in a simi- 
lar way as B Xg'y The only difference is 



that A„ for 6 — > ^7 is not small relative to A* and 
Ac; therefore, also the current-current operators 
O" and O2, weighted by A„, contribute. Un- 
fortunately, these operators induce long-distance 
contributions to B ^ X^j, which at present only 
can be estimated using models. In ref. |3^, these 
long-distance effects were absorbed into the the- 
oretical error. 

Using fib = 2.5 GeV and the central values 
of the input parameters, the analysis in ref. |3^ 
yields a difference between the LO and NLO pre- 
dictions for BR{B — > Xd'y) of ~ 10%, increas- 
ing the branching ratio in the NLO case. For a 
fixed value of the CKM-Wolfenstein parameters p 
and rj, the theoretical uncertainty of the average 
branching ratio {BR{B Xdj)) of the decay 
B Xd'y and its charge conjugate B Xdj is: 
A{BRiB ^ Xd"f))/{BR{B ^ ^^7)) - ±(6 - 
10)%. Of particular theoretical interest for con- 
straining p and rj is the ratio of the branching 
ratios, defined as 



Rld'y/s^) 



^ {BR{B ^ Xd-f)) 
~ {BR{B ^ Xgj)) ' 



(2.11) 



in which a good part of the theoretical uncer- 
tainties cancels. Varying the CKM-Wolfenstein 
parameters p and 77 in the range —0.1 < p < 0.4 
and 0.2 < rj < 0.46 and taking into account 
other parametric dependences, the results (with- 
out electroweak corrections) are 

6.0 X 10"*^ < BR{B Xdi) < 2.6 x 10"^ , 
0.017 < Rid-f/sj) < 0.074 . 

Another observable, which is also sensitive to the 
CKM parameters p and 77, is the CP rate asym- 
metry acp, defined as 



acp 



T{B ^ Xd-f) - T{B ^ Xdj) 
T{B ^ Xdj) + T(B ^ X~d-f) 



(2.12) 



Varying p and 77 in the range specified above, 
we obtained 7% < acp < 35% We would 
like to point out that ac p is at the moment only 
available to LO precision and therefore suffers 
from a relatively large renormalization scale de- 
pendence. 

A measurement of the semileptonic FCNC 
decays B Xgl+l- and B Xdi+i^, below 
the J/t/"- and above the p, w-resonance regions 
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in the dilepton invariant mass, can also be used 
to extract |Vfs| and \Vtd\, respectively. In this 
context, these decays and the related ones, B 
Xsvv and B — > Xdvv, were discussed some time 
ago iQ. The decays B {Xs,Xd)vv are prac- 
tically free of long-distance contributions |^ and 
the renormalization scale dependence of these de- 
cay rates has also been brought under control 
||3^ . Hence, these decays are theoretically re- 
markably clean but, unfortunately, they are dif- 
ficult to measure. The ALEPH collaboration has 
searched for the decay B —t Xgi^D, setting an up- 
per boimd BR{B X^vv) < 7.7 x 10"-* (at 90% 
C.L.) which is a factor 20 away from the SM 
expectations 

In contrast, the prediction of the decay rate 
T{B Xsi^d") still suffers from many uncer- 
tainties. The most important ones are due to 
intermediate cc states. Because of the non per- 
turbative nature of these states, the differential 
spectrum can be only roughly estimated when 
the invariant mass mi+i- is not sufficiently be- 



low 



However, for low s = m^^^^_ jm 



relatively precise determination of the spectrum 
is possible using perturbative methods only, up 
to calculable HQET corrections. The dominant 
HQET corrections were evaluated ^ and found to 
be smaller than 6% for 0.05 < s < 0.25. There- 
fore, the B Xsi~i^ rate integrated over this 
region of s should be perturbatively predictable 
as precisely as T{B Xg'y), i.e. up to about 
10% uncertainty. However, the presently avail- 
able NLO QCD corrections |||, ^ have not yet 
reached this precision. The formally LO term 
is suppressed, which makes it as small as some 
of the NLO contributions. Consequently, some 
NNLO terms still can be large. The NNLO pro- 
gram for this decay was recently started |4^, by 
calculating the two-loop matching conditions for 
all the operators relevant for B Xgl^t^ ■ The 
improved matching allows to remove an impor- 
tant ±16%) matching scale uncertainty for 
BR{B Xst^l~) in the mentioned region of 
s, leading to BR{B -> Xs£+i~) = 1.46 x lO^^. 
A remaining perturbative uncertainty of about 



^For an overview of different treatments of the 
and l/m^ corrections in i? — > {Xs,X^)i^ , and of the 
^^Qjy/'i^ terms, generated by the «-quark loops in _B ^ 
Xal+l-, we refer to 



13%, due to the unknown two-loop matrix ele- 
ments of the four-quark operators, was also esti- 
mated in ref. lIlJ. 



3. 5 ^ in generalized two-Higgs 
Doublet models 

Two Higgs Doublet Models (2HDMs) are con- 
ceptually among the simplest extensions of the 
SM. Studies of BR{B X^^i) in these models 
can already test whether the observed high ac- 
curacy of the NLO SM result is a generic feature 
of NLO calculations [|[ |2^ or a rather particluar 
one, valid for the SM only. Such studies can obvi- 
ously provide also important indirect bounds on 
the new parameters contained in these models. 

The well-known Type I and Type II mod- 
els are particular examples of 2HDMs, in which 
the same or the two different Higgs fields cou- 
ple to up- and down-type quarks. The second 
one is especially important since it has the same 
couplings of the charged Higgs to fermions 
that are present in the Minimal Supersymmetric 
Standard Model (MSSM). The couphngs of the 
neutral Higgs to fermions have important differ- 

li. How- 



ences from those of the MSSM |44 



ever, since beside the W, only charged Higgs 
bosons mediate the decay B Xg^ when ad- 
ditional Higgs doublets are present, the predic- 
tions of BR{B Xs-i) in a 2HDM of Type II 
give, at times, a good approximation of the value 
of this branching ratio in some supersymmetric 
models |||. 

It is implicit in our previous statements that 
we do not consider scenarios with tree-level flavour 
changing couplings to neutral Higgs bosons. We 
do, however, generalize our class of models to 
accommodate Multi-Higgs Doublet models, pro- 
vided only one charged Higgs boson remains light 
enough to be relevant for the process B Xg^- 
This generalization allows a simultaneous study 
of different models, including Type I and Type II, 
by a continuous variation of the (generally com- 
plex) charged Higgs couplings to fermions. It al- 
lows also a more complete investigation of the 
question whether the measurement of BR{B 
Xg'y) closes the possibility of a relatively light 

not embedded in a supersymmetric model. 
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We will show that the NLO QCD corrections 
to the Higgs contributions to BR{B Xsj) are 
much larger than the corresponding corrections 
to the SM contribution ||^ , irrespectively of the 
value of the charged Higgs couplings to fermions. 
This feature remains undetected in Type II mod- 
els, where the SM contribution to BR{B X^^) 
is always larger than, and in phase with, the 
Higgs contributions. In this case, a comparison 
between theoretical and experimental results for 
BR{B — > Xs^) allows to conclude that values of 
m^-i = 0{mw) can be excluded. Such values 
are, however, still allowed in other 2HDMs. 

These issues are illustrated in Sec. 3.3, after 
defining in Sec. 3.1 the class of 2HDMs consid- 
ered, and presenting the NLO corrections at the 
amplitude level in Sec. 3.2. 

3.1 Couplings of Higgs bosons to fermions 

Models with n Higgs doublets have generically a 
Yukawa Lagrangian (for the quarks) of the form: 

-K '^1 4, - K h U'j,^ + h.c. , (3.1) 



where g^, </>,;, [i — 1, 2) are SU(2) doublets {(jii — 
«(T^(/)*); w'^, d'^ are SU(2) singlets and h'^, /i" 
denote 3x3 Yukawa matrices. To avoid flavour 
changing neutral couplings at the tree-level, it is 
sufficient to impose that no more than one Higgs 
doublet couples to the same right-handed field, 
as in eq. (3.1). 



After a rotation of the quark fields from the 
current eigenstate to the mass eigenstate basis, 
and an analogous rotation of the charged Higgs 
fields through a unitary n x n matrix U, we as- 
sume that only one of the n—1 charged physical 
Higgs bosons is light enough to lead to observable 
effects in low energy processes. The rt-Higgs dou- 
blet model then reduces to a generalized 2HDM, 
with the following Yukawa interaction for this 
charged physical Higgs boson denoted by H^ : 



g ( Y- M 
—= <XulV 

V2 L mw 



YuR 



nil, 



■Vdr\H 



mw 



(3.2) 

In (3.2), V is the Cabibbo-Kobayashi-Maskawa 
matrix and the symbols X and Y are defined 
in terms of elements of the matrix U (see cita- 
tions in ref. j|]). Notice that X and Y are in 



general complex numbers and therefore poten- 
tial sources of CP violating effects. The ordinary 
Type I and Type II 2HDMs (with n = 2), are spe- 
cial cases of this generalized class, with {X, Y) = 
(— cot /3, cot /3) and {X,Y) ~ (tan /3, cot /3), re- 
spectively. 

3.2 NLO corrections at the amplitude level 

It turns out that the charged Higgs contributions 
do not induce operators in addition to those in 
the SM Hamiltonian i?]^^ in eq. (p.ll). More 



specifically, only step 1) below eq. ( |2.2D , gets 
modified when adding the charged Higgs boson 
contributions to the SM one. The new contri- 
butions to the matching conditions have been 
worked out independently by several groups ||47| , 
^ ^ , by simultaneously integrating out all heavy 
particles, W, t, and iJ+ at the scale fiw- This is 
a reasonable approximation provided 'mfj± is of 
the same order of magnitude as mw or m^. 

Indeed, the lower limit on m//± coming from 
LEP I, of 45 GeV, guarantees already = 
Olmw)- There exists a higher lower bound from 
LEP II of 55 GeV for any value of tan/3 g for 
Type I and Type II models, which has been re- 
cently criticized in ref. This criticism is 
based on the fact that there is no lower bound 
on iriA and/or uih coming from LEP |9|. 

After performing steps 1), 2), and 3) listed 
below eq. (^), it is easy to obtain the quark 
level amplitude A{b 57). As the matrix ele- 
ments {s^\Oi\b) are proportional to the tree-level 
matrix element of the operator O7, the amplitude 
A can be written in the compact form 



A{b S7) 



4Gf 
V2 



VtbV,:D{sj\07\b)t 



(3.3) 

For the following discussion it is useful to decom- 
pose the reduced amplitude D in such a way that 
the dependence on the couplings X and Y (see 
eq. (|3.2|)) becomes manifest: 



D = DsM+XY*DxY 



\Y\^Dyy 



(3.4) 



In Fig. ^ the individual D quantities are shown 
in LO (dashed) and NLO (sohd) order, for mH± = 
100 GeV, as a function of fib', all the other input 
parameters are taken at their central values, as 
specified in ref. |^]. To explain the situation, one 
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Figure 2: LO (dashed) and NLO (solid) predictions 
of the various pieces of the reduced amplitude D for 
m^± = 100 GeV (see text). 

can concentrate on the curves for Dxy ■ Starting 
from the LO curve (dashed), the final NLO pre- 
diction is due to the change of the Wilson coeffi- 
cient C7, shown by the dotted curve, and by the 
inclusion of the virtual QCD corrections to the 
matrix elements. This results into a further shift 
from the dotted curve to the solid curve. Both 
effects contribute with the same sign and with 
similar magnitude, as it can be seen in Fig. ^. 
The size of the NLO corrections in the term Dxy 
1) is 



AD 



■NLO 



XY 



D 



■LO 
XY 



D 



LO 
XY 



D 



LO 
XY 



-40%! (3.5) 



A similarly large correction is also obtained for 
Dyy- For the SM contribution Dsm, the situa- 
tion is different: the corrections to the Wilson co- 
efficient C7 and the corrections due to the virtual 
corrections in the matrix elements are smaller in- 
dividually, and furthermore tend to cancel when 
combined, as shown in Fig. |^ 

The size of the corrections in D strongly de- 



pends on the couplings X and Y (see eq. (|3.4|)): 
AD/D is small, if the SM dominates, but it can 
reach values such as —50% or even worse, if the 
SM and the charged Higgs contributions are sim- 
ilar in size but opposite in sign. 

3.3 Results 

The branching ratio BR{B ^s7) can be sche- 
matically written as 

BR{B X,-f) cx \D\^ + ■■■ , (3.6) 



where the ellipses stand for Bremsstrahlung con- 
tributions, electroweak corrections and nonper- 
turbative effects. As required by perturbation 



theory, \D\'^ in eq. (| 



\D\ 



\D 



LO, 



^) should be understood as 
AD 



1 + 2Re 



-LO, 



-LO 



(3.7) 



i.e., the term \AD/D^^\'^ is omitted. If Re(AD/ 
D ) is larger than 50% in magnitude and nega- 
tive, the NLO branching ratio becomes negative, 
i.e. the truncation of the perturbative series at 
the NLO level is not adequate for the correspond- 
ing couplings X and Y. As shown in ref. j^, this 
can happen even for modest values of X and Y. 

However, theoretical predictions for the branch- 
ing ratio in Type II models stand, in general, on 
a rather solid ground. Fig. ^ shows the low-scale 
dependence of BR{B Xg'y) for matching scale 
^iw — 'rnfj± , for 'mfj± > 100 GeV. It is less than 
±10% for any value of m^± above the LEP I 
lower bound of 45 GeV. Such a small scale un- 




100 200 300 400 500 600 

Figure 3: BR{B Xs-y) in a Type II model with 
tan /3 = 2, for various values of Hb. The leading QED 
corrections are included (see text). 

certainty is a generic feature of Type II models 
and remains true for values of tan /3 as small as 
0.5. In this, as in the following figures where 
reliable NLO predictions are presented, the re- 
cently calculated leading QED corrections are in- 
cluded in the way discussed in the addendum ||25|] 
of ref. 1^. They are not contained in the result 
shown in Fig. ^, which has an illustrative aim 
only 

In Type II models, the theoretical estimate of 
BR{B — > Xs^) can be well above the experimen- 
tal upper bound of 4.5 x lO""* ( 95% C.L.) |7|, 
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leading to constraints in the (tan (3, mH± ) plane. 
The region excluded by the CLEO bound, as well 
as by other hypothetical experimental bounds, is 
given in Fig. ^. For tan/3 = 0.5,1,5, we exclude 



800 
600 
400 

200 








a.ono-i 










: \ 






3.5*10-* J 




4.0*10-* 






" 1 


4.5*10-* ~_ 

... 1 



10 



100 



tan /S 



Figure 4: Contour plot in (tan f3, mu± ) in Type II 
models, obtained by using the NLO expression for 
BR{B Xsj) and possible experimental upper 
bounds. The leading QED corrections are in- 
cluded. The allowed region is above the correspond- 
ing curves. 



respectively m^± < 280, 200, 170 GeV, using 
the present upper bound from CLEO. 

Also in the case of complex couplings, the re- 
sults for BR{B Xsj) range from ill-defined, 
to uncertain, up to reliable. One particularly in- 
teresting case in which the perturbative expan- 
sion can be safely truncated at the NLO level, 
is identified by: Y = 1/2, X = 2exp(i(/)), and 
mH± — 100 GeV. The corresponding branching 
ratios, shown in Fig. ^ are consistent with the 
CLEO measurement, even for a relatively small 
value of mH± in a large range of (j). Such a light 
charged Higgs can contribute to the decays of the 
t-quark, through the mode t — s- H^b. 

The imaginary parts in the X and Y cou- 
plings induce -together with the absorptive parts 
of the NLO loop- functions- CP rate asymmetries 
in _B ^ Xs^. A priori, these can be expected to 
be large. We find, however, that choices of the 
couplings X and Y which render the branching 
ratio stable, induce in general small asymmetries, 
not much larger than the modest value of 1% ob- 
tained in the SM Is^. 



x" 

A 

ICQ 
K 

m 



_ 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 _ 




NLO - 






1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



50 



100 



150 



Figure 5: BR{B Xs"/) as a function of (f), 
where (p parametrizes X — 2exp(i<^), for Y — 1/2, 
m.ff± — 100 GeV. Solid, dashed and dash-dotted 
lines correspond to = 4.8, 2.4, 9.6 GeV. The lead- 
ing QED corrections are included (see text). Super- 
imposed is the range of values allowed by the CLEO 
measurement. 

4. X,7 in SUSY models 

Rare decays also provide guidelines for super- 
symmetry model building. Their observation, 
or the upper limits set on them, yields strin- 
gent constraints on the many parameters in the 
soft supersymmetry-breaking terms. The pro- 
cesses involving transitions between first and sec- 
ond generation quarks, namely FCNC processes 
in the K system, are considered to be most ef- 
ficient in shaping viable supersymmetric fiavour 
models. 

The severe experimental constraints on flavour 
violations have no direct explanation in the struc- 
ture of the MSSM. This is the essence of the well- 
known supersymmetric flavour problem. There 
exist several supersymmetric models (within the 
MSSM) with speciflc solutions to this problem. 
Most popular are the ones in which the dynam- 
ics of flavour sets in above the supersymmetry 
breaking scale and the flavour problem is killed 
by the mechanisms of communicating supersym- 
metry breaking to the experimentally accessible 
sector: In the constrained minimal supersym- 
metric standard model (mSUGRA), supergrav- 
ity is the mediator between the supersymmetry 
breaking and the visible sector ||5^ . In gauge me- 
diated supersymmetry breaking models (GMSBs), 
the communication between the two sectors is re- 
alized by gauge interactions |5|]. More recently. 
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the anomaly mediated supersymmetry breaking 
models (AMSBs) were proposed, in which the 
two sectors are linked by interactions suppressed 
by the Planck mass |5^. Furthermore, there are 
other classes of models in which the flavour prob- 
lem is solved by particular flavour symmetries. 

Neutral flavour transitions involving the b 
quark, do not pose yet serious threats to these 
models. Nevertheless, the decay B Xg^ has 
already carved out some regions in the space of 
free parameters of most of the models in the 
above classes (see and references therein). 
In particular, it dangerously constrains several 
somewhat tuned realizations of these models . 
Once the experimental precision is increased, this 
decay will undoubtedly help selecting the viable 
regions of the parameter space in the above class 
of models and/or discriminate among these or 
other possible models. It is therefore important 
to calculate the rate of this decay as precisely as 
possible, for generic supersymmetric models. 

As we saw in section 2, BR{B Xs^) is 
known up to NLO precision in the SM. The cal- 
culation of this branching ratio within general 
supersymmetric models is still far from this level 
of sophistication. There are several contributions 
to the decay amplitude: Besides the SM- and 
the charged Higgs one, there are also chargino- 
, gluino- and neutralino contributions. All these 
were calculated in within the mSUGRA model. 
The inclusion of LO QCD corrections was as- 
sumed to follow the SM pattern. A calculation 
taking into account solely the gluino contribution 
has been performed in Js^ for a generic super- 
symmetric model, but no QCD corrections were 
included. 

An interesting NLO analysis of B ^ Xs^f 
was recently performed Q in a specific class of 
models where the only source of flavour violation 
at the electroweak scale is encoded in the CKM 
matrix. The calculations were done in the limit 

/ig = 0{mg,mq,mij :^ nw 

= 0{mw,mH±,m^±,m^^), (4.1) 

and terms of order [fiw I l^gY ip > 2) were dis- 
carded. At the scale /ivi/ the new contributions 
can be matched onto the same operators as in the 
SM (see also ref. @). It is shown in the analysis 



Q that much lower values for m[{± are allowed 
than in the typc-II 2HDM, discussed in section 2, 
due to the possiblity of destructive interference 
between the charged Higgs and the chargino con- 
tributions. It is illustrated in Fig. ^ that in such a 
cancellation scenario the NLO QCD corrections 
are important: The uppermost curve is the LO 
result for the type-II 2IIDM with mH± = 100 
GeV; switching on also the chargino contribu- 
tion at LO (using the parameters mentioned in 
the caption) , leads to the second curve (from the 
bottom). The NLO result for the combination 



20 



LOcomplete 
LOdecoupled 
NLO - 
2HDM-L0 
2HDM-NL0 



0.4 



0.6 0.8 1 
heavy masses (TeV) 



Figure 6: BR(B —> Xs^) as a function of fig for the 
parameters: tan/3 — 1, mjf± = = m^^ ~ 100 



GeV, rriy 



300 GeV, Or 



-n/lO, At 



At; 



all heavy particle masses equal to ^g; the lighter 
chargino is predominantly higgsino. Figure taken 
from M. Ciuchini et al. ||| 

of the charged-Higgs and the chargino contribu- 
tion is represented by the middle curve. 

This calculation, however, cannot be used in 
particular directions of the parameter space of 
the above listed models in which quantum ef- 
fects induce a gluino contribution as large as the 
chargino or the SM contributions. Nor it can 
be used as a model-discriminator tool, able to 
constrain the potentially large sources of flavour 
violation typical of generic supersymmetric mod- 
els. 

The flavour non diagonal vertex gluino-quark- 
squark induced by the flavour violating scalar 
mass term and trilinear terms is particularly in- 
teresting. This is generically assumed to induce 
the dominant contribution to quark flavour tran- 
sitions, as this vertex is weighted by the strong 
coupling constant g^. Therefore, it is often taken 
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as the only contribution to these transitions and 
in particular to the B —^ Xs"f decay, when ex- 
tracting order of magnitude upper bounds on 
flavour violating terms in the scalar potential 
jS^ , [STt - Once the constraints coming from the 
experimental measurements are imposed, how- 
ever, the gluino contribution is reduced to values 
such that the SM and the other supersymmet- 
ric contributions cannot be neglected anymore. 
Any LO and NLO calculation of the B Xgj 
rate in generic supersymmetric models, there- 
fore, should then include all possible contribu- 
tions. 

The gluino contribution, however, presents 
some peculiar features related to the implementa- 
tion of the QCD corrections. In ref. |^ this con- 
tribution to the decay b ~^ sj is therefore investi- 
gated in great detail for supersymmetric models 
with generic soft terms. It is shown that the rela- 
vant operator basis of the SM effective Hamilto- 
nian gets enlarged to contain magnetic and chro- 
momagnetic operators with an extra factor of as 
and weighted by a quark mass mi, or mc, and also 
magnetic and chromomagnetic operators of lower 
dimensionality, as well as additional scalar and 
tensorial four-quark operators. A few results of 
the analysis in ref. j^] are given in the following, 
showing the effect of the LO QCD corrections on 
constraints on supersymmetric sources of flavour 
violation. 

To understand the sources of flavour viola- 
tion which may be present in supersymmetric 
models in addition to those enclosed in the CKM 
matrix, one has to consider the contributions to 
the squark mass matrices 



supersymmetry breaking potential. As a conse- 
quence, gluino contributions to the decay 6 — > 57 
are induced by the off-diagonal elements of the 
soft terms mjj^^, m) ,^^., ^),lr and 

It is convenient to select one possible source 
of flavour violation in the squark sector at a time 
and assume that all the remaining ones are van- 
ishing. Following ref. |5^, all diagonal entries in 
"^d,LL' "^^RW and rn\ j^i^ are set to be equal 
and their common value is denoted by m?. The 
branching ratio can then be studied as a function 
of 



TO? 



RR,ij 



TO? 



'JLR,ij 



("^ LR)ij 



(4.4) 
(4.5) 



The remaining crucial parameter needed to de- 
termine the branching ratio is a; = to? /to?, where 
Trig is the gluino mass. In the following, we con- 
centrate on the LO QCD corrections to the gluino 
contribution. In Figs. and |[ the solid 
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M} 



'^f,RL ""^fMR 



(4.2) 



(4.3) 



Ff,LL + Df^LL Ff,LR. 

Ff,RL Ff,RR + ^f.RR . 

where / stands for up- or down-type squarks. In 
the super CKM basis where the quark mass ma- 
trices are diagonal and the squarks are rotated in 
parallel to their superpartners, the F terms from 
the superpotential and the D terms turn out to 
be diagonal 3x3 submatrices of the 6x6 mass 
matrices Al^. This is in general not true for the 



additional terms (4^), originating from the soft 



Figure 7: Gluirio-iriduced branching ratio BR{B 
Xs'y) as a function of a; = m|/m?, obtained when the 
only source of flavour violation is 5lr,23 (see text). 

lines show the QCD corrected branching ratio, 
when only Slr,23 or 5ll.23 are non vanishing. 
The branching ratio is plotted as a function of x, 
using TOg = 500 GeV. The dotted lines show the 
range of variation of the branching ratio, when 
the renormalization scale Hh varies in the inter- 
val 2.4-9.6 GeV. Numerically, the scale uncer- 
taintly in BR{B X^^) is about ±25%. An 
extraction of bounds on the 5 quantities more 
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Figure 8: Same as in Fig. |^ when only 5ll,23 is 
non- vanishing. 



precise than just an order of magnitude, there- 
fore, would require the inclusion of NLO QCD 
corrections. It should be noticed, however, that 
the inclusion of the LO QCD corrections has al- 
ready removed the large ambiguity on the value 
to be assigned to the factor as(/^) in the gluino- 
induced operators. Before adding QCD correc- 
tions, the scale in this factor can assume all val- 
ues from 0{mb) to 0{m\Y)' the difference be- 
tween BR{B — > Xs^) obtained when as{mb) or 
when as{mw) is used, is of the same order as the 
LO QCD corrections. The corresponding values 
for BR{B Xg'^) for the two extreme choices 
of ji are indicated in Figs. ^ and ^ by the dot- 
dashed lines (/i = mw) and the dashed lines 
(/i — rrib). The choice /i = mw gives values 
for the non-QCD corrected BR{B — > Xsj) rel- 
atively close to the band obtained when the LO 
QCD corrections are included, if only 23 is 
non-vanishing. Finding a corresponding value of 
fj, LhaL mliiimlzes Lhe QCD correcLlous in Llie case 
studied in Fig. |^, when only Slr,23 is different 
from zero, depends strongly on the value of x. 
In the context of the full LO result, it is impor- 
tant to stress that the explicit as factor has to 
be evaluated - like the Wilson coefficients - at a 



5. Summary 

Significant progress in the theoretical description 
of rare B decays has been achieved during the 
last few years. NLO QCD corrections are avail- 
able for radiative inclusive decays in the SM. 
Power correction (1/m^, 1/m^) are also under 
control. The description of the Fermi motion of 
the b quark in the B meson has been refined. Im- 
portant NNLO QCD improvement of the match- 
ing conditions of the operators relevant for B — > 
Xs£~^£~ was made, which removes the ±16% match- 
ing scale uncertainty in the invariant mass dis- 
tribution of the lepton pair. In some regions of 
the parameter space, NLO QCD corrections to 
BR{B Xsj) in 2HDMS are huge. They are, 
however, under control in the typc-II model and 
are important to derive reliable bounds on tan (3 
and mjf± . NLO QCD corrections to BR{B 
Xg'y) are available also in particular SUSY sce- 
narios. Calculations for more general situations 
are in progress. 
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